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Abstract 

We outline a microscopic approach to the superconducting fluctuations and 
pairing correlations in the attractive Hubbard model above T c , using the func- 
tional integral method. A crossover from BCS superconductivity to Bose 
condensation of preformed pairs is studied by constructing the appropriate 
Ginzburg-Landau functionals. The differences between the lattice and the 
continuum models are discussed. The case of quasi-two-dimensional super- 
conductors as well as the model with non-local pairing interaction are also 
examined. The effects of gaussian fluctuations of the order parameter are 
analyzed in the T-matrix approach, the self-consistent Hartree approach and 
finally within the Ginzburg-Landau theory. Simlarities with a paramagnon 

theory of itinerant-electron magnetism are pointed out. 
PACS 74.20.-z 
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I. INTRODUCTION. 



The problem of a crossover from a weak coupling BCS superconductivity to Bose con- 
densation of preformed local pairs has recently attracted great attention, stimulated by the 
experimental results cocerning unusual properties of new high T c superconductors as well 
as the discovery of Bose condensates in optically trapped dilute atomic gases. It has been 
established that high temperature superconductors (HTS) generally exhibit a low carrier 
density, small value of the Fermi energy (oc 0.1 — 0.3 eV), short coherence length £ and 
they are extreme type II superconductors with the Ginzburg-Landau parameter k ^> 1 [1-4]. 
Experimental studies of HTS and also other non conventional superconductors indicate that 
the cuprate high T c , bismuthates, fullerenes and Chevrel phases belong to a unique group 
of superconductors characterized by high transition temperatures relative to the values of 
n s /m* . These materials have their T c proportional to Tp (the Fermi temperature) or Tb (the 
Bose- Einstein condensation temperature) [2], with T B ~ (3 — 30)T C and T F pz (10 — 100)T C . 
Moreover, the existence of the pseudogap phase in the underdoped HTS is well experimen- 
tally established and it can be interpreted as a state with pairing without long-range phase 
coherence [4,5]. The above points show that these materials are not in a weak coupling BCS 
limit, but rather in the intermediate to extreme strong coupling regime. 

One of the central questions regarding the physics of HTS is the understanding of the 
evolution from a weak coupling limit of large Cooper pairs to the strong coupling regime 
of small local pairs with increasing coupling constant and description of the intermediate 
(crossover) regime. 

Several attempts to describe the crossover have employed the BCS mean field theory, 
generalized to allow a change in the chemical potential due to the pairing away from the 
weak coupling limit [6-11]. It has been shown that the BCS wave function for T = OK is a 
reasonable ansatz for any attraction strength. For the fermions in continuum, the controlling 
parameter is £,okp (£o - the coherence length , kp - the Fermi momentum) and the BCS 
and BC limits correspond to £o^f 1 and £o^f <C 1, respectively. An extension to finite 
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temperatures has been first proposed by Nozieres and Schmitt-Rink [9] , who included single- 
pair fluctuations within a T-matrix approximation and found that T c for the 3D continuum 
dilute fermion system evolves smoothly between the limits (see also Refs. [12-14]). Similar 
analysis has also been carried out for a 2D dilute attractive Fermi gas [15,16]. 

In this paper we overview the functional integral approach to the crossover problem 
using mainly the Hubbard model with on-site attraction. The attractive Hubbard model 
is a generic model of superconductivity on a lattice which can describe the evolution from 
the BCS type superconductivity to a Bose condensation of real space pairs, with increasing 
coupling strength (for review see Refs. [1,17,18]). The layout of this paper is the following. 
In Section II, we present application of the functional integration method to the negative 
U Hubbard model. In Section III, the crossover is studied by a detailed derivation of the 
Ginzburg-Landau theory in the weak attraction limit and the corresponding Ginzburg- 
Pitaevskii functional in the strong coupling regime. We discuss the differences between 
the lattice and continuum fermion models with s-wave attraction. The case of quasi-2D- 
superconductivity and the model with pair hopping interaction are also discussed. In Section 
IV, we discuss the Gaussian fluctuations and the T-matrix approach. The self-consistent 
approach to the superconducting fluctuations and pairing correlations is further discussed 
in Section V, where a derivation of a self-consistent Gaussian theory is given. In Section VI 
we analyze the amplitude fluctuations above T c using the Ginzburg-Landau theory. 



II. FUNCTIONAL INTEGRAL APPROACH. 

Here we will outline the functional integral approach to the attractive Hubbard model. 
The Hamiltonian of the model can be split into two parts 

H = H + H 1} (1) 

#0 = Yl tiA C 3° ( 2 ) 

#i = -|£/|£pjp;, (3) 



where p\ = cLcji , and Uj = Uj — /iSij. Uj is the transfer integral, |C/|-the on-site attraction 
and \i stands for the chemical potential. The grand partition function of the model (1) in 
the interaction representation with respect to H is given by 

Z = Trexp[-P(H + HJ] = Tr[e- pH °T T exp(- [ drH^r))} = 

Jo 

Z (r T eM- f Q 'drH^r))^ , (4) 

where T T is the "time" ordering operator, i?i(r) = exp(rH )Hi exp(— tH ),Z = 
Tre~P Ho , (•••)o = Tr(---e-P H »)/(Tre~P Ho ) and (3 = (hsT)' 1 . The functional integral 
representation of Z in terms of auxiliary fields Aj(r) is obtained in a standard manner by 
making use of the Hubbard-Stratonovich transformation [19] 



Z = Z j D 2 Aexp(-^[Ai(r)]), 



(5) 



where 



fi[A] = - In ( T T 



*[A i (r)] = 5:/" /, dr|A i (r)| 2 /|^| - n[A], 
i J0 

exp(Y, j*dT[k i (T)p\(j)+H.c]) 



(6) 
(7) 



D 2 A = DAD A*. Introducing the Nambu representation = (ctpCjjJ, one obtains for 

n[A] 



ft [A] = - In (T T 



(8) 



where A,- 



A* 



The matrix Green's function 



^■(ry) = -(T T * i (r)*}(^)> 



f G tJ (r,T>) F tJ (r,T>) ^ 
4(r,r') -G tJ (r,T>) 



\ ij 



(9) 
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where the field operators are given in the Heisenberg representation with respect to the 
Hamiltonian (1), can be evaluated within the functional integral scheme as follows 



g i Ar,r') = (g t ,(T,r',{A})) A , (10) 

where the average (• • -) A = / _D 2 Aexp(— \P[A])(- ■ ■)/ / D 2 A exp(— *[A]) denotes the func- 
tional average. The Green's function: Qij(r, r', {A}) is just the Green function for electrons 
in the random fields Aj(r) and is given as 

e w (r,r,{A})= ^-y- 

= -(T T Mr)m) H , +Hr ■ ("I 
S(A)=T rCT p('-Ajf' 3 dr*tA,(r)*,(T)j , (12) 

(• • -) Ho+Hr denotes the average with the Hamiltonian H + H r with H given by Eq.(l) and 
H r given below 

H r = -\J2(plMr) + H.c.) , (13) 

i 

and A is a formal parameter being set to 1 at the end. This local Green's function satisfies 
the BCS equation in the presence of spatially and time varying local order parameter: 

E (~ - *>* " aA ^)) Gij(t, t', {A}) = lSijSir - t') , (14) 

where a z is the Pauli matrix. An exact representation for the partition function can be given 
by 

Z = Z Q J D 2 Aexp[-V] 
* = Y1 f^dr\A t (r)\ 2 /\U\ - Tr ln(l - AQ ), (15) 



where 



G 1 = (--^ - Uja^J 5ij5(T - t') 



(16) 



We note that formal power series expansion of \& will contain only even powers in A and it 
has the form 

c^Tr (AQo) 

Trln(l-Ag ) = -E \ k ' ■ (17) 



The explicit form of the expansion up to the 4th order is given below 



* (4) =E(l^ + x(q^n)) |A q (^ 

q,n \ \ U I / 



+^7J3 E r+ + (Pi,P2,P3;n 1 ,n 2 ,n 3 )A* pi (ie ni )A* p2 (ie n2 ) 

A p3 (l€ n3 ) A pi -|_p 2 _p 3 (l£ ni + is n2 — i£ n3 ) , (18) 

In the above equation 

is the free particle pairing susceptibility, /(x) = (exp(f3x) + l)" 1 is the Fermi distribution 
function, e n = jn are the bosonic Matsubara frequencies. A^(ie n ) is the Fourier transform 
of the field Aj(r) and is given by 

f/3 



A q (^n) = ^ E / Q dre^-^JA^r). (20) 
The quartic term is of the form 



T + + (Pi, P2, Ps]n 1 ,n 2 ,n 3 ) = J— ^g^{k,iu m )g {k + p 1 ,iuj m -ie n ) 

Go (k + p 3 - p 2 ,iuj m + ie n2 - ie ns )G ~(k + pz,iw m - ie ns ), (21) 

where Qq 1 = iu n — (e k — fi)o z is the Fourier transform of the free fermion Green's function 
and uj n = jf(n + \) are the fermionic Matsubara frequencies. If we neglect all the frequency 
dependence of A in the expansion(17) (the static approximation), then we have the following 
result for the quartic term (we also neglect the momentum dependence): 

2 

, e k = e k - fj, , (22) 
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(iu m ) 2 ~ £k 2 _ 

which is the standard expression for the quartic term in the Ginzburg-Landau (G-L) theory. 
We note that the corresponding coefficients can be evaluated explicitly and basing on such 
expansion one can study the question of the crossover from the BCS pairing (when fi Ep) 
to the bosonic limit of tightly bound pairs (fx ps —E b /2, where E b is the binding energy for 
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two fermions in an empty lattice). For completeness, we give below the time-independent 
Ginzburg-Landau functional 



q V / fci,fc2,fc3 



&3 J 



(23) 



where A q = A q (0). The coefficients appropriate to a weak coupling superconductor are 
given by 

n-N(F)\n T r - J^Ei^l h - 7 C(3)iV(^) 
a-N(E F )\n-,c- , b - , (24) 

where N(E F ) is the density of states (DOS) at the Fermi level and E F is the Fermi energy. 
The above G-L functional will be used in Section VI in discussing the effects of amplitude 
fluctuations above T c . 

III. FROM BCS TO BOSE SUPERCONDUCTIVITY 

As mentioned above, given the functional integral representation (15), the evolution from 
BCS to Bose type superconductivity can be analyzed. The coefficient in quadratic term of 
(18) is of the form 

A(q, ie n ) = |^ + x(q, i£ n)- ( 25 ) 

It is convenient to eliminate the pairing potential \U\ in favor of the two-particle binding 
energy in vacuum, E b , [7], through the relation 

A—E^— (26) 
\U\ N^2e k + E b 1 ; 

In the long-wave limit, assuming the parabolic spectrum for fermions with the mass m, the 
expansion is given by 

q 2 

A(q, ie n ) = a + c- ide n , (27) 

2m 

where 
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1 y 1 1 y> tanh(/3e fc /2) 

a N^2e k + E b N 2 ^ 2e k ' 1 } 

1 ^ r Aanh(/3e fc /2) /3 \ 



W T V 4 ^ 8£ fe cosh 2 (/3£ fe /2); 
/^tanh^£fc/2) /? _ 2 tanh(/3e fc /2) \ (k ■ n) 2 

+ V. cosh 2 (/^ fc /2) + e fe cosh 2 (^V2) 4 / 8me* 1 J 

6/2 l V Aanh(fe/2) g 1 \ 

7 *VV 44 84cosh 2 (^/2)J' ( j 

and the coefficient J will be specified later. 

In the BCS/Fermi liquid regime E b / Ep <C 1, ~ £f ^ 7" and we recover the standard 

form of the Ginzburg Landau functional: 



(31) 



q = (q,ie n ), with the coefficients a,b,c given by Eqs.(24) and with d = nN(E F )/8kBT c . 
Indeed, with the rescaling $(&) = \/2cA k (the Gorkov relation), it takes the conventional 
form of the time-dependent G-L functional. 

In the opposite Bose limit E b /Ep ^> 1, and the chemical potential— \x approaches E b /2 
at the transition point, thus —pL ^> T c . Defining the effective chemical potential // = p — 
we then obtain, a = -^2p,c = ^,d= fg 1 , where N(0) is the DOS in 2D. After the 
rescaling of the pair-field amplitude $(q, ie n ) = \j ^^-A n (ie n ), the effective action becomes 

Seff = E - 2/i - <en) l<%)| 2 + 6 E (2)$(3)$(1 + 2-3), (32) 

Q \ / 1,2,3 

where m* = 2m, p = E b /2 + /i, and b is the parameter describing the interaction between the 
tightly bound pairs. Eq.(32) yields required Gross-Pitaevski functional describing interacting 
bosons in the dilute limit, with the effective mass m* = 2m and the effective chemical 
potential p [20,21]. The parameter b describes repulsive interaction between the bosons and 
it strongly depends on dimensionality. For a constant DOS, the above procedure yields 

V2 - (33, 
lb 

where N(0) = In the limit -> E b /2, the rescaling gives b m 2/N(0). Thus, the 
only interaction is due to the Pauli principle for constituent fermions. However, in general, 



one would expect, that the effective boson-boson repulsion should also be controlled by the 
scattering length for two bosons. In the strictly d=2 case, the scattering amplitude at long- 
waves and low-frequencies vanishes, so the above mapping must be taken with care, and 
the effects of higher orders in the expansion (18) may be of importance. We also notice 
that in 2D, the coefficient b of the GL functional (Eq.(33)), when ji — > 0, is divergent, thus 
indicating a singular point in the evolution from BCS to Bose limit. This however can be 
again cured by taking into account the interaction between the fluctuations [22]. 

A question of the derivation of the time-dependent Ginzburg-Landau functional by this 
method, in the weak and strong coupling regimes, is discussed in Ref. [12]. 

Finally, we should mention that the effective mass of the bosons m* in the continuum limit 
approaches 2m. If the tight binding spectrum of the form e k = —2t(cos(k x a) + cos{k y a) + 
cos(k z a)), a being the lattice constant, is used, one obtains the following expansion for 
A(q, ie n ): 

A(q, ie n ) = a + Cq 2 — ide n) (34) 
where a and b are given by Eq.(28) and Eq.(31) and 

f n are the derivatives of the Fermi function. The coefficient d can be obtained by first 
making an analytical continuation %(q = 0, ie n ) — > x(q = 0, u + i0 + ) and then expanding 
in uj [12]. One gets: 

- 1 _ 1 - 2f(e k ) . ttAW x ( ^ 

where 6(/i) is the Heaviside step function. In the Bose limit the coefficients take the form: 
a = —2p,/U 2 ,C = 2t 2 a 2 /\U\ 3 ,d = 1/U 2 , and upon rescaling the pair-field amplitudes one 
gets the functional (32), with the bosonic mass m* = \U\/ (4t 2 a 2 ). This result is in accord 
with the perturbation expansion for the attractive Hubbard model in the strong coupling 
limit [1]. 
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A. Quasi-two dimensional case 



Let us now consider the quasi-two dimensional case and assume the following spectrum 
for the fermions interacting via the on-site attractive potential U : 

k 2 



2m i 



— 2t± cos(k z d) : 



(37) 



where k 2 = k 2 + k 2 , m\\ is the in - plane fermion mass and t± is the hopping amplitude 
between the neighboring planes at distance d. In contrast to the strictly 2D case, now the 
s-wave BCS pairing and the bound state formation are not necessarily related to each other. 
The bound state condition Eq.(26), shows that for a given ratio of t\\/t± the threshold must 
be reached for \U\, to form the bound state in vacuum. In the BCS case, evaluation of 
the quadratic and quartic coefficients will yield an appropriate generalization of the G-L 
functional to the quasi 2D case. 

Let us consider the strong coupling case above the threshold to form the bound pair, and 
where \x is below the bottom of the fermionic band and the temperature region ksT < E b . 
In such a case we can neglect the Fermi factors in the evaluation of the quadratic coefficient 
and obtain 



1 1 
A(q, is n ) = - £ £ _ At± cos(M) + At± + Eb 



N T ^ + £l - 4*± cos(q z d/2) cos(k z d) + 4t ± - 2/i - ie n ' 

and q = (q,q z ). 

After integrating over k, one gets: 



(38) 



A(q,ie n ) = ^^ln 



4mu 



V 



2//-i£ n + y(£ ? -2//-i £n ) -16ticos 2 (^) 
E b + 4t ± + ^{E b + 4t ± ) 2 - ml 



(39) 



/ 



Introducing the effective chemical potential by the relation \i = —\E\> — 2t± + fi, and 
expanding (39) to lowest order in q\\ and ie n , we obtain 



10 



A(q, ie n ) = ^ (■£- + ^(1 - cosferf)) -2H- ie n ) . (40) 



By rescaling the pairing-field amplitude $(q, ie n ) = y ^^-A q (i£„), one obtains the 
Gross-Pitaevski action for the quasi-2D-bosons 

= E + mW 1 - cos( ** d » - 2p " fe ») l$M|2 + 

6^$'(l)$'(2)$(3)$(l + 2-3), (41) 

1,2,3 

where m_i_ = l/t5_<i 2 is the fermionic mass in the perpendicular direction. The bosons have 
the in-plane mass 2m\\ and the perpendicular mass m^j_ = 2E b m 2 ± d 2 being much greater 
than m±. The physical origin of the increased perpendicular mass is clear, to move the pair 
between the planes one has to break it first. Indeed, these (composite) bosons will condense 
at the temperature given by [23] 

= ™*d „ Ef (4) 

Bc 2m^\n{2k b T c vE b m\d i ) ~ 2\n{{Ak b T c E b m\d A Y V ; 

n* = n/2 is the boson density, n being the electron density and Ep = nnd/m\\. Thus, 

the transition temperature for the quasi-two dimensional case in the preformed pair regime 

is proportional to Ep. This temperature will decrease with increasing interaction making 

the pair binding stronger. We should point out that independently, same result has been 

obtained in Ref.( [24]), by using a different approach. 

B. The extended Hubbard model with pair-hopping interaction 

The above derivations of the Ginzburg-Landau functional can be extended to the case 
of the Hubbard model with the pair hopping interaction (the Penson-Kolb-Hubbard model), 
which involves the non-local pairing interaction and is given by the Hamiltonian: 

H = E Uj^Cja + P \ Pl - X - YX J aP\Pi + H - c -) > ( 43 ) 

i,j,<r i i,j 

where Jy is the pair- hopping (intersite charge exchange) parameter [25]. The functional 
integral representation of the partition function of the model (43) is given by Eq.(15) with 
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* = E |Ag( ' £w)|2 - Trln(l - A&), (44) 

q,n 

where Vq = —U + J q and J q = 2J(cos(k x a) + cos^a) + cos(A; z a)) is the Fourier transform 
of Jij. The coefficient of the quadratic term in the power series expansion of ^ in A is 

A(q, ie n ) = tt + x(q, *£n), (45) 

where the pairing susceptibility is determined from Eq.(19). The small q, u expansion is 
now of the following form 

A(q, ie n ) = ai + Ciq 2 - ide n , (46) 

where 

1 1 v tanh(/3e fc /2) 1 _ 1 1 _ tanh(/3e fc /2) 

Ci = % + C, (47) 

and C is given by Eq.(35). The complete expansion of the functional (44) to the 4th order 
is fully specified by the above expressions together with the coefficient J given by Eq.(36) 
and b given by Eq.(30). 

IV. GAUSSIAN FLUCTUATIONS AND THE T MATRIX APPROACH 

In this Section we will analyze the contribution from the Gaussian fluctuations to the 
free energy and the self-consistent T-matrix approach. The effect of Gaussian fluctuations of 
the order parameter is obtained by keeping the quadratic terms only in the full functional. 
The corresponding contribution to the thermodynamic potential is simply given by 

/^ (2) = ^£ln[l + |^|x(q,^n)], (48) 



qn 



where 



_1 

x(q, ie n ) = -~tz £ G (k, iu m )G (q - k, ie n - iu m ), (49) 

• km 
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with G (k, iu n )~ l = iu n — {e k — /i), and is given by Eq.(19). 

The BCS pairing instability corresponds to a singularity for q = 0, uj = in JF^ 2 ), and 
the condition (the Thouless criterion) 1 + |C/|x(0, 0) = gives T c . 

Equations (48,49) can be used to form a fully renormalized, conserving approximation 
for the free energy and Green function [26-28]. Let the functional be defined as 

$[G] = ^£ln[l-|^lx(q^n)], (50) 
X(q, ie n ) = G ( k > iu m )G(q - k, ie n - iu m ), (51) 

" km 

where G(k, ioj n ) is the full Green function and we use conventional definition of the pairing 
susceptibility Eq.(51). 

The corresponding self-consistent approximation for the free energy is given by the func- 
tional 

2 

•FpjG] = -—J2exp(iuj n r]) {£(k, iu n )G(k, iu n ) + ln[-G (k, iio^ 1 + E(k,iu; n )]} (52) 

• kn 

evaluated at its stationary point with respect to variations of G and E. At this stationary 
point G, E, and $ are related by: 

G(k, iu n ) = [G {k, iuJnY 1 - E(k, iuj n )Y l , (53) 

(54) 

Evaluation of E(k, iuj n ) yields 

£(k, iu n ) = -j— T(q, ie m )G(q - k, ie m - iu n ), (55) 

qm 



where 



T(q,ie n ) = - — J^l , (56) 

i - l^lx(q,^„) 



is the T-matrix. 
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The electron concentration as a function of T and /i is obtained as n(T,fi) = —dJ-'/dfi. 
By using stationary properties of T one gets the expected result 

2 

n(T, fi) = -Q^^2 exp(iu n r))G(k, iu n ). (57) 

" kn 

The above equations (53,55,56,51,57) constitute the system of self-consistent T-matrix equa- 
tions and they were recently analyzed numerically for the attractive Hubbard model in two 
dimensions (mostly in low concentration regime) [29]. 

V. SELF-CONSISTENT APPROACH TO THE SUPERCONDUCTING 

FLUCTUATIONS 

A. Hartree theory and variational principle approach 

We would like to derive the best Gaussian form for the superconducting fluctuations. 
Let us seek for a trial functional ^> T in the following form 

*T[A q ]=j:A q \A q \ 2 , (58) 
i 

where q = (q, ie n ). Application of the Feynman variational principle allows us to write the 
corresponding variational free energy F Q 

(3F = - In J D 2 Ae~ w ^ ] + (* - ^ T ) T , (59) 

where (• • -) T = J D 2 A exp(-* T [A,])(- • •)/ / D 2 A exp(-* T [A,]). To evaluate the RHS of 
Eq.(59) we consider two cases: 

1. We use expansion of the exact functional ^ truncated at the 4th order term. 

2. We use the exact form of * given by Eq.(15) 

In the first case we use \1/ as given by Eq.(18) and after exploring the properties of Gaussian 
integrals we obtain for F 
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pF = - In y D 2 Ae~* T ^ + £ 



— + - A q 



2iV/3 



Er + - + -( P1 , P2 , P2 )(|A P1 | 2 )(|A P2 | 



Pl,P2 



£r + - + -( Pl ,p 2 , Pl )<|A P1 | 2 ><|A P2 | 



(60) 



2^ P t P2 

where (|A p | 2 ) is the functional average with The best coefficients A q are obtained from 



the stationary condition <5 ^°^ 



0. With the notation t q =< \A q \ 2 > for the mean square 



amplitude of the order parameter (fluctuating in space and time), we obtain 



A q = + x(q) + E [r + ~ + ~(^,P) + r + ~ + -(p, g,p) 



tp. 



(61) 



Since \I/t is the Gaussian functional then £ 9 = and we obtain the following self-consistent 
equation for t q : 



(62) 



If we neglect the frequency dependence and perform small momentum expansion for x(q) 
and also neglect the momentum dependence of the T terms in (62), then we recover the 
Hartree equation for the amplitude fluctuations used in Section VI. 

B. Fully self-consistent approach 



In the second case we use the complete functional (15) to evaluate the trial free energy 
(59) [31]. By making the substitution A q = A q l l 2 ^ q , one gets 



(63) 



The stationary condition 
for A a 



5JJ3Fq) 



0, after integrating by parts, yields the following expression 



A 



5 2 ^ 



(64) 
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i.e. the optimal quadratic functional coefficients are self-consistently averaged second deriva- 
tives of the true functional. We note that the random phase approximation (RPA) result is 
obtained if the second derivative is evaluated at the origin: 

5 2 ^ 



SA q 5A* q 



A=0, 



(65) 



and it simply yields Eq.(25) with the free particle-particle susceptibility (19). By using 
Eq.(64) and the free energy functional Eq.(15), after tedious manipulations (similar to those 
in Ref. [30]), one arrives at 

—^(-TrMl-AG,)) = £^'{A}^, fe , +9 {A}, (66) 

1 1 k,k' 

where Q^,{A}, a = 1,2, denotes the diagonal matrix element of the space-time Fourier 
transform of ftj(r,r',{A}) given by Eq.(14). We note that G{A} = 0„(1 - A^ ) _1 is 
the electron propagator in the random potential Aj(r). The expression for the coefficient 
A q = A(q,ie n ) is then given by 

A(q,ie n ) = ^ + x(q,ie n ), (67) 

where 

*(q,ien) = --^ E (^ 11 (k,k / ,^ m ;{A})^ 22 (k + q,k / + q,^ m+ „;{A})) T (68) 

^ k,k',m 

Therefore, we have derived a self-consistent theory of superconducting fluctuations in which 
the paring susceptibility x(k, (ie n ) is given by the averaged particle-particle propagator in 
the presence of the random fields. In this approach the particle-particle propagator in 
Eq.(67) is just the product of two local Green functions self-consistently averaged with the 
Gaussian weight exp(— A q \ A q \ 2 ). (To lowest order in \U\, when we ignore the random fields, 
X reduces to the bare Xi an d we recover the RPA result). The above theory is very similar in 
the spirit to the paramagnetic polaron theory in itinerant-electron magnetism [30] and can 
be used to describe incoherent local pair state. Finally, we note that the presented theory 
when resorting to the static approximation yields a theoretical foundation for the approach 
in which the electronic DOS is evaluated with the assumption of existence the local order 
parameter amplitude above T c and averaged with the Gaussian weight. 
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VI. AMPLITUDE FLUCTUATIONS ABOVE T c IN THE GINZBURG-LANDAU 
THEORY AND INCOHERENT LOCAL PAIR STATE 

It is of interest to study further the fluctuation effects (and pairing correlations) and 
their role in establishing pair formation above T c , within the conventional Ginzburg-Landau 
theory. From the previous analysis, one comes to the conclusion that above T c there will 
be a region of parameters in which there are pairs with amplitude and phase fluctuating 
in space and time. Such a state, however, has no long-range phase coherence which is 
necessary to establish the superconducting state. This scenario is somewhat similar to a 
local-moment itinerant electron magnetism in transition metals and alloys, in which the 
magnetic moments exist above transition temperature but their direction is random. A 
possibility of such a state in the negative — U Hubbard model has been considered earlier 
[1] to describe the intermediate region between BCS and the local pair regime. Such an 
incoherent local pair state has also been analyzed by Gyorffy et al [32,33] within the coherent 
potential approximation to the negative —U Hubbard model with the use of a constant 
density of states. 

In order to show the existence of a finite gap amplitude above T c , we consider the 
standard Ginzburg-Landau functional (G-L) and apply the Hartree approximation to treat 
the interaction between superconducting fluctuations [33]. 
In zero magnetic field one has 



F[*(r)] = j 



dr 



aft)|^(r)| 2 + &|^(r)| 4 + |^|W(r)| 2 



(69) 



where t = T/T c and m is an effective electron mass. The superconducting order parame- 
ter ^(r) can be related to the gap parameter A(r), for example, via the Gorkov relation: 
tf(r) = ( 

8(IisT c ) 2 )^ 2 ^( r )' w ith n the number of electrons per unit volume. With the Fourier 
transform to reciprocal space we have 



F[A k ]/fi = £ 



[a{t) + cA; 2 ]|A k | 2 + b £ A k A p A q A*_ r 



p+q 

k L p,q 



(70) 
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where Q is the volume of the system. Assuming that b and c are temperature independent 
and performing the simplest quadratic factorization of the last term we obtain the gaussian 
form of the G-L functional: 



^[AJ/n = E 



k L 



a(t) + ck 2 + 2bJ2 < \ A q\ 2 > 



lAwl 2 . 



(71) 



In the above equation the expectation value has to be calculated as the functional average: 

/D 2 A(...)exp[-F[A k ]/(A; B T)] 



(....) = 



(72) 



JD 2 Aexp[-F[A k ]/(k B T)] " 

The mean square gap amplitude A 2 = < |Ak| 2 > is then given by the following self 
consistent equation 



k B T 



E 



i 



a(t) + ck 2 + 2bA 2 



In 2D we obtain 



(73) 



ln(l + w/x) = (d/t)(x + l-t), 



where the following relations have been used 



a{t) = a(0)(l -t),t = T/T c , |A | 2 = -a(0)/26, 

x = A 2 /|A | 2 -l + t,w = (k c £ ) 2 , 
J _ 4 7 rla(0)lg|A | 2 7C(3)ft 2 n| A | 2 



k B T c 



87rm(k B T c 



In 3D case we get 



w 

w* - V^arctan(^) = (d*/t)(x + 1 - t), 



where 



w* = k c £ , d* = 



7((3)h 2 n\A \% 



(74) 



(75) 



16m(k B T c y 

£o and A are the zero temperature coherence length and the zero temperature energy gap, 
respectively, k c denotes the momentum cutoff. 
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We have solved these equations with the use of parameters appropriate for YBa 2 Cu307, 
taking |A | = 20meV, T c = 92 K, n = 6 * 10 14 cm~ 2 (for 2D case) and n ~ 10 21 cm~ 3 (for 
3D case), m = 5m e , k c = ^q 1 . The result is illustrated in Figs. 1-2 for several values of d 
and d* . For both 2D and 3D cases for T above T c , A smoothly decreases and saturates at 
a finite value. Depending on the values of d and d*, A is of order(0.1 — 0.2)A for 2D case 
and of order (0.05 - 0.125)A in 3D. 

The parameter k c is not precisely defined in the Ginzburg-Landau theory and we consider 
it to vary. If we allow randomness of A(r) down to the atomic scale, then k c can be taken 
as an inverse of the lattice constant (4A), giving w = 9 (2D) and w* = 3 (3D). In this case, 
we observe (Figs. 1-2) that after an initial decrease, A/A starts to increase and saturates at 
larger values. The gap magnitude is of the order (0.15 — 0.275)A in 2D and (0.1 — 0.2)A 
in 3D, depending on the values of the parameters d and d*. Thus if we lower k c , increasing 
temperature is more effective in making larger thermal fluctuations in the amplitude of 
the local order parameter A(r). This characteristic behavior is very similar to that known 
in the theory of local moment magnetism [34]. The above calculations carried out at the 
phenomenological Ginzburg-Landau level clearly illustrate the importance of the fluctuations 
above T c and point out toward a possibility of existence of a state with incoherent pairs. In 
the incoherent local pair state the one-particle density of states will develop a pseudogap or 
real gap depending on the electron density and the strength of the attractive interaction. 
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FIG. 1. versus T/T c in two dimensions for two different values of the cutoff k c . 
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FIG. 2. ^ versus T/T c in three dimensions for two different values of the cutoff k c 
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